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We present calculations of the upper critical field H C 2(T) in superconducting films as a function 
of increasing disorder (as measured by the normal state resistance per square Ra). In contradiction 
to previous work, we find that there is no anomalous low temperature positive curvature in H C 2 
vs T as disorder is increased. We show that the previous prediction of this effect is due to an un- 
justified analytical approximation of sums occuring in the perturbative calculation. Our treatment 
includes both a careful analysis of first-order perturbation theory, and a non-perturbative resumma- 
tion technique. No anomalous curvature is found in either case. We present our results in graphical 
form. 



I. INTRODUCTION 

Increasing-disorder is known to suppress superconductivity in low-dimensional systems such as thin films and 
narrow wiresaii This occurs because the disorder causes electrons to move diffusively rather than ballistically, making 
them less efficient at screening the Coulomb repulsion between electrons. The increased Coulomb repulsion decreases 
both the electron-electron attraction needed for superconductivity, and the density of states of electrons available 
for pairing at the Fermi surfaceO. Typical types of (experimental data are: (i) T C (R^), the transition temperature 
as a function of normal state resistance per squareaiMj; (ii) Ao(-Rq), the order parameter at zero temperature, as a 
function of normal state resistance per squareOEj; (in\-.Mr&CT , Ra) = T c (Ra,H), the upper critical field as a function 
of temperature and normal state resistance per squarefltfrEj; (iv) T c (-Rc, 1/t s ), transition temperature as a function 
of resistance per square and spin-flip scattering rate in films with magnetic impurities.EJ It is found experimentally 
that T c (i?n) curves from a wide variety of materials fit a universal curve with a single fittiftg,|piLrameter, whilst the few 
experimental measurements of Ao(i?o), seem to have Ao(R a )/2kBT c (Rc) roug hly constantly. This fitting to a single 
curve, whilst pleasing in showing that the basic ingredients of our theories are correct, does not allow detailed analysis 
of the theory. Data of types (iii) and (iv) are more promising because there is an additional parameter to vary - the 
magnetic field in (iii), and spin- flip scattering rate in (vA- To the best of our knowledge, only one experiment of. tY i P f 
(iv) has been performedMl, and we discuss it elsewhere. t£l Several experiments of type (iii) have been performedBE-3E-3; 
some seem to show a positive curvature in H C 2(T) at low temperature as disorder is increased. Moreover this effect 
is predicted by theoryOlla, and this seems to be another confirmation of the basic theoretical model. However, we 
need to be careful: positive curvature in H C 2(T) is a ubiquitous feature of exotic superconductorsEj, and occurs in 
many systems where localization is not believed to be the cause. Indeed any pair-breaking mechanism that varies as 
a function of magnetic field can lead to such anomalous behaviour in H C 2- This means that it is often difficult to 
distinguish between the various mechanisms that might be present. It is therefore particularly important to be sure 
of our theory, and in this light we re-examine the predictions of localization theory. 

One of the main problems of the localization theory is that even first-order perturbation theory results are hard 
to obtain correctly. The first-order results are capable of explaining experimental data in the weak disorder regime, 
but for strcuiger disorder it is clear that we need something else. As an example consider the prediction for T c 
suppression!™, 

where T c o is transition temperature for clean system, Rq — 2iTh/e 2 ss 162/cf2, and r is the elastic scattering time. This 
yields an exponential curve for T c (i?n), which behaves like a straight line for small R^. It is clear that T c deduced from 
this equation can never go to zero for finite Ru, as happens in experiment. A very simple ad hoc way of going beyond 
simple perturbation theory is to replace T c o on the right hand side by T c , pleading perhaps to self-consistency. If we 
define x — In (T c0 /T c ), (3 = In (1/2ttT c0 t), and t = R^/Rq, the new equation has the cubic form 

x= t -{fi + xf (2) 

and can easily be solved. However a new problem emerges because there are two positive roots for every value of Ru. 
At first we can take the larger of the roots, on physical principles, because it is this root which tends to T c /T c0 = 1 



f 



at Ra — 0. However we eventually come to a re-entrance point beyond which no solutions exist. It is clear that this 
re-entrance is unphysical, an artefact of our ad hoc extension of perturbation theory. In the case of T c (Ru) the story 
has a happy ending in that perturbation theory can be correctly extended by a renormalization group (RG) treatment 
based on Finkelstein's interacting non-linear sigma modeloEd. This leads to the result 

J!!)'- >J11*\ (3) 
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where 7 = — 1/ln (T c o/1.13r). This equation reduces to the first-order result for small t and now T c goes smoothly to 
zero at t = 7-. The three curves are plotted for comparison in Fig. (1). 
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FIG. 1. Theoretical predictions for transition temperature as a function of resistance per square. The solid line is the 
first-order prediction; the dashed line is the "self-consistent" first-order prediction, which suffers from re-entrance; the 
long-dashed line is the result obtained from either renormalization group or ladder re-summation approaches. 

The reason for discussing TJJia) in detail above is that the same problem occurs for H C 2(T,Ra) = T c (R a ,H). 
The standard theory in this caseO, due to Maekawa, Ebisawa and Fukuyama (MEF) is the equivalent of the ad hoc 
extension discussed above, and has the form 
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where 1/th = DeH. This equation suffers similar re-entrance problems at finite 1/tjj to those found at 1/th = 0, 
when it is just the T c equation. Indeed the H C 2(T) curves can only be plotted down to the value of T at which 
re-entrance occurs, and at this point the curves appear to have infinite slope. This leads to us asking the question 
of whether the positive curvature in H C 2(T) is also an artefact of the ad hoc approximation used. What we need 
to answer this question is the finite magnetic field analogue of the RG result discussed above. However the RG is 
very difficult, and the answer is not forthcoming from this source. Fortunately Oreg and FinkelsteinEiJ have recently 
shown that the RG result can be obtained from a diagrammatic resummation technique. This method has the great 
advantage of being analytically tractable and easy to use. In this paper we extend this approach to finite magnetic 
field to see whether we really do expect positive curvature in H C 2(T). Indeed this paper is intended as somewhat of 
a showcase for this resummation method, to demonstrate the ease of its extension to a wide variety of problems. 

Oreg and Finkel'steinEd take the Coulomb interaction to be featureless, and we shall follow them. We should 
therefore explain why it is legitimate to use a featureless interaction rather than the correct screened Coulomb 
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interaction. The screened Coulomb interaction has a singularity at low momentum, and one might naively think that 
this would lead to a strong enhancement of the suppression of transition temperature. However a cancellation-occurs 
between diagrams 1-4 and diagram 5 of Fig. 2, which effectively removes the low momentum singularitytra. The 
Coulomb interaction is then effectively featureless, all diagrammatic sums being dominated by large frequency and 
momentum. It turns out that one gets the same result from doing the perturbation theory correctly with the screened 
Coulomb interaction as from using a constant interaction of strength g = N(0)V = 1/2 and keeping only diagrams. 
1-4. This is what we shall do in this paper. We also note that the re-summation method of Oreg and FinkelstcinEil 
keeps only diagrams 3 and 4, which is again legitimate as their contribution is greater than that from diagrams 1 
and 2. However it is not difficult to include these terms in the re-summation, as we will show in section III. The 
cancellation of low momentum singularities is expected to be a general feature, enforced by gauge invariance, so we 
are free to ignore them and use a featureless interaction as long as we include the diagrams that give the dominant 
contribution at large frequency and momentum. 
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FIG. 2. The diagrams needed to calculate the first-order perturbative correction to the pair-propagator L(q, 0). (a) Definition 
of pair propagator L(q, 0) in terms of the pair polarization bubble P(q, 0) and the 4-point BCS interaction — \>y\. (b) Zeroth-order 
(mean-field) polarization bubble in a dirty superconductor. The solid lines are electron Green functions; the dashed lines are 
impurity lines, (c) First-order perturbative correction to pair polarization bubble. The wiggly line is the screened Coulomb 
interaction. 
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The outline of the rest of the paper is as follows. In section II we review the first-order perturbation calculation for 
H C 2 (T) assuming a featureless Coulomb interaction. We derive the analytical MEF formula by making an asymptotic 
approximation to the Matsubara frequency sums which arise in perturbation theory. We make plots of H C 2 (T) both by 
solving the implicit MEF equation, and its equivalent where the Matsubara sums are performed exactly. Surprisingly 
we find no positive curvature in the latter case. In section III we derive H C 2(T) using the resummation technique 
both in the form used by Oreg and Finkel'stein, and also in an extended form which includes self-energy diagrams 
missed in their formalism. Again we find no positive curvature in H e 2(T). In section IV we discuss the experimental 
situation and draw conclusions. 



II. REVIEW OF FIRST-ORDER PERTURBATION THEORY 

In this section we will carefully review the first-order perturbation theory calculation^] of H C 2 (T) . We do this in 
detail because we will show that accurately performing sums rather than making an analytical approximation to them 
removes the upward curvature in H C 2(T). We identify T c as the temperature at which the magnetic field dependent 
pair propagator diverges. To find the latter we calculate the correction to the momentum-dependent pair propagator, 
L(q,0), and then make the usual substitution Dq 2 — > 2/th, where 1/th = DeH. The zeroth-order (mean field) pair 
propagator is given by 



LoMr 1 = N(0) 
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so that upon substitution Dq 2 — » 2/tjj we get the usual Abrikosov-Gorkov resultS for the pair-breaking effect of the 
magnetic field 



v T c a J V 2 / \ 2 2itT c t h , 
We will calculate the corrections to the pair polarization bubble, SP(q,0), which will lead to a change in T c given by 

Since we will assume a featureless interaction, we need to (evaluate diagrams 1 to 4 of Fig. 2. Diagrams 1 and 3 
involve the summation of three terms to form a Hikami boxcJ, the general form for which is 

2^V(0)r 4 [-D(A\ + 2Al) + |ei| + |e 2 | + |e 3 | + |e 4 |] 6(-e 1 e 2 )6(-e 2 e ;i )6(-e 3 e i ) (8) 

4 2 

where A4 = J]] q« • qi+i and A 2 = q; • q;+2, the being the incoming momenta, and the e, the Matsubara 

i=l i=l 

frequencies on the electron Green functions in the box. Using standard diagrammatic rules, and the above result for 
the Hikami box, we obtain the 4 contributions to the pair polarization bubble 

p a MrtvrSrrSrS^ Dg 2 + Dg' 2 + 3^1 + \e n \ 
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We note that the relative signs of e; and e„ are irrelevant in the sum, since the summand depends only upon |ej| and 
|e n |, and that the featurelessness of the potential allows us to shift the momentum sum in terms 2 and 4. We find 
that terms 1 and 2 partially cancel, whilst terms 3 and 4 reinforce, to yield the sum 
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where g = N(0)V and we have performed the g'-sum and set Dq 2 — > 2/th- 

We first reproduce MEF's analytic approximation to the sums over Matsubara frequencies. To do this it turns 



out to be easier to make a choice of relative sign of Matsubara frequencies, eje„ < 0, and to set e„ = e; 
ei = 2irT(l + 1/2) and uj m — 2TrTm, the sum becomes 
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where a = 1/2-kTth, M — \/2ttTt. If we first evaluate the sum for a — 0, we see that it will be dominated by large 
to, so we replace the difference of digamma functions by In to, and the sum over to by an integral to get 
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To evaluate the result for finite a we subtract off the a = result and try to analytically approximate the difference. 
Again we expect the sum to be dominated by large to, so we can ignore the a in ij)(l/2 + m + a) and (to + 2a), and 
ignore the term proportional to a as it is less divergent. The difference then has the form 
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From our knowledge of the calculation which includes the full screened Coulomb interaction we know that we should 
set 3 = 1/2, from which it follows that g/4ir 2 N(0)D = t. Putting the Eqns. @ and @ into Eqn. (0) yields the 
MEF formula cited in Eqn. (g). It turns out that this approximation, although apparently reasonable, is not justified. 
If we plot the values of P(a) calculated by performing the sums directly, we find that they do not agree with Eqn. (Q). 
Consequently the H C 2(T) curves predicted by the analytic approximation and exact sum are also different. 




FIG. 3. Plots of H C 2(T) deduced from first-order perturbation theory with "self-consistency". The curves on the left are 
derived using the analytic MEF formula, and clearly show upward curvature at low temperature as resistance is increased. The 
curves on the right are derived by calculating Matsubara sums exactly, and show no upward curvature. It follows that the 
upward curvature is an artefact of making an incorrect analytical approximation in the derivation of the MEF formula. Both 
curves terminate at non-zero values of temperature because of re-entrance problems. 
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On the left side of Fig. 3 we display H C 2 (T) predicted by the MEF result of Eqn. (Q) for several values of i? n 
in the case where \ti(1/2t:T c qt) — 6. We clearly see the positive curvature at low temperature, and the result that 
the curves all terminate at finite T due to re-entrance problems. On the right side of Fig. 3 we plot H C 2(T) deduced 
directly from the first-order perturbation theory result of Eqn. (|ll|). We note that since the answers we get depend 
upon the logarithm of the upper cut-off, if we treat the upper cut-off differently, we will get slightly different answers. 
However, this difference will always be lost in fitting to experiment since we determine the unknown upper cut-off 
parameter from the initial slope of the data. The reason we mention this is that the exact sums can be performed in 
two slightly different ways: we can have the sum over ei and u> m , or the sum over q and e„, the cut-offs in each case 
being taken at 1/r. The results are slightly different due to different treatment of the upper cut-off. However in both 
cases we find no positive curvature in H C 2(T). There is still the problem that the curves finish at finite T because 
of re-entrance, and this tells us that the ad hoc prescription being used to go beyond first-order perturbation theory 
is still inadequate. But again we stress the key point of this section: even within the ad hoc extension of first-order 
perturbation theory, we do not get positive curvature in H C 2 (T) if we do the Matsubara sums exactly. 



III. OREG AND FINKELSTEIN'S RESUMMATION TECHNIQUE 

In this section we answer the question of what is the correct way to go beyond first-ordex. perturbation theory. 
To do this, we extend a resummation technique recently developed by Oreg and FinkelstcinEd (OF). This involves 
calculating the pair scattering amplitude, r c (e n ,e;), and identifying T c as the temperature at which it first diverges. 
The ladder summation involved is shown diagrammatically in Fig. 4 This is actually an extended version of the 
OF approach because as well as considering diagrams 3 and 4 of Fig. 2, which correspond to an effective Coulomb 
pseudopotential, and end up in the block tA, we also consider diagrams I and 2 of Fig. 2, which correspond to an 
effective self-energy, and end up in the block £ which renormalises the Cooperon impurity ladder. In fact we will 
consider both versions of the summation technique: the simpler one involves using the bare Cooperon Co rather 
than the dressed one C. OF demonstrate that the matrix equation they obtain in 2D can be approximated by a 
parquet-like differential equation that turns out to be identical to that obtained from renormalization group analysis. 
Instead of generating a continuous approximation to the matrix equation, we simply solve the equation and change 
the temperature (which affects the upper cut-off M = 1/2itTt) until we find the first zero eigenvalue of the matrix, 
at which point it has become singular, and we have reached T c . Numerically this involves diagonalizing matrices of 
rank less than 2000 or so. 
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FIG. 4. Diagrammatic equation for the scattering amplitude matrix T(e„, e;). Block 7 is the BCS interaction. Block tA is the 
correction to the effective interaction caused by the interplay of Coulomb interaction and disorder. Block E is the correction 
to the bare Cooperon ladder Co caused by Coulomb interaction and disorder. We consider both the original resummation 
procedure which uses the bare Cooperon Co, and an extended procedure which includes the fully dressed Cooperon C. 



Let us now proceed to the calculation of the pair amplitude matrix T. If we evaluate the diagrams of Fig. 4, but 
ignore the magnetic field and self-energy correction to the Cooperon, we obtain the same equation as OF, 
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except that we have explicitly kept both positive and negative Matsubara frequencies because the expression for 
A(q, e„) will turn out to depend upon whether its two Matsubara frequencies have the same or opposite sign. This is 
a reflection of the breaking of time-reversal invariance caused by the magnetic field. Let us first see how to solve the 
above equation, and then later put in the magnetic field, and the correction to the Cooperon. As a matrix equation, 
it takes the form 



f = -| 7 |i + iA-i[-| 7 |i + tA]r 1 f 

where T nm = T(e„, e TO ), l nrn = 1, A nm = A(e„, e m ) and e nm = (n + l/2)5 nm . This has the solution 
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and I nm = 5 nm is the identity matrix. It follows that when the matrix II has an eigenvalue equal to I/I7I, the pair 
amplitude diverges, and we have found T c . Our approach is therefore to start at the mean field transition temperature, 
T c o, and decrease temperature until one of the eigenvalues of (J— |7|II) changes sign. The matrix II (T) depends upon 
T both through the dependence of A upon T and through its rank 2M, where M = (2ttTt)^ 1 . We start at the 
mean-field value of M, which we will call Mq, and decrease the temperature by increasing M successively by one. We 
diagonalize the matrix II for each value of M until an eigenvalue changes sign. At this point we have found T c for 
the given problem, and T c /T c q — Mq/M. We can then change a parameter such as t or a and repeat the procedure. 
This removes the need for any perturbative expansion. 

Let us now put back in the magnetic field and Cooperon self-energy corrections. The |e m | denominator in the 
Eqn. (14) comes from the Cooperon 
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and in the presence of a magnetic field, \e m \ is replaced by \e rn \ + 1/th, and hence e nrn = (n + 1/2 + a)5 nm , where 
a — (2itTth The contributions to tA from the diagrams of Fig. 4 which correspond to diagrams 3 and 4 of Fig. 2 
are given by 
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and so performing the g'-sum and setting Dq 2 
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It follows that the matrix elements of A are 

In 



M 



hi 



71 + 771 + 1 

M 



2a 



71 + 771 + 1 J (n + 771 + 1) 



e n e m > 



(21) 



To include the self-energy correction into the Cooperon, we note that the corrected Cooperon is given by 
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so that we need to absorb a factor — 1/2ttN(0)t 2 into S. The contributions from the diagrams of Fig. 4 which 
correspond to diagrams 1 and 2 of Fig. 2 are 
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which partially cancel to give the result 
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The weak-localization contribution to the Coopcron self-energy is given by 
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Incorporating £ into the Cooperon means that e has elements 
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Using the new formulas for e and A, we now plot H c2 (T) for various values of t = Bv/Rq. The results are shown 
in Fig. 5: the plot on the left does not include self-energy corrections; the plot on the right does. The results are 
very similar and show no sign of upward curvature in H c2 (T) . We note that there are no re-entrance problems in our 
calculation of T c : we can plot the curves down to as low a temperature as we like if we are prepared to diagonalize 
large enough matrices. 
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FIG. 5. Plots of H C 2(T) deduced from the non-perturbative resummation method. The curve on the right includes the 
self-energy corrections to the Cooperon ladder; the curve on the left does not. Neither curve show any sign of positive 
curvature at low temperature. 



IV. DISCUSSION AND CONCLUSIONS 



The central message of the paper is that the theory of localization and interaction does not predict an anomalous 
positive curvature in the upper critical field H c2 {T) of thin film superconductors. A subsidiary message is that the 
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resummation method developed by Oreg and FinkelsteinEll is a very powerful and adaptable tool for going beyond 
perturbation theory in disordered superconductors. We suspected that the positive curvature in H C 2(T) was an artefact 
related to the re-entrance problem found in T c (Ra), and that it would not survive a systematic non-perturbative 
treatment. The latter is true, but the cause of the positive curvature artefact turns out to be an incorrect analytic 
approximation to sums of Matsubara frequencies. 

We will now consider the experimental situation by attempting to fit experimental data to both the MEF formula 
of Eqn. , and the exact first-order perturbation theory result of Eqn. ([ll]) . The data which fits best to the localization 
theory is that of Graybeal and Beasleyu on amorphous films of Mo-Ge, because values of T c for all films can be obtained 
by choosing a single value of the parameter (3 = In (1/2ttT c0 t). The fits in Fig. 6 used the parameter value (3 = 7.17 
for the MEF curves and (3 = 5.3 for the exact sum curves. Values of T c q — 7.2K and H c q — 120kOe were taken 
directly from the data. The experimental data appears to fit better to the exact sum curves than the MEF curves, 
and there seems little sign of an upward curvature in H C 2(T). 

When the same procedure is applied to the data of Okuma et alli3 on Zn films, we find that we need to use 
different values of (3 for the two films to get the correct T c for given R^. The plots in Fig. 7 use values (3 = 7.5 for the 
400ft film and (3 = 8.1 for the 600ft film for the MEF formula; (3 = 5.56 /or the 400ft film and (3 = 6.1 for the 600ft film 
for the exact sum formula. Both films have the same thickness of 100A, but different resistances, and hence different 
diffusion constants D. The latter were calculated using the material parameters for Zn to get values of H C 2(T = 0) 
for zero-resistance films of the same composition equal to 5.08kOe for the 400ft film and 7.62kOe for 600ft film. T c o 
was taken to be IK. Again the data seems to-fit better to the exact sum curves than the MEF formula. 

To fit the data of Hebard and Paalanenta on amorphous In-InO^ films again requires a different value of (3 for 
each film. The plots in Fig. 8 use (3 = 3.5 for the 2250ft film and (3 = 4.0 for the 2900ft and 3300ft films for the 
MEF curves; f3 = 1.6 for 2250ft film and (3 = 2.0 for the 2900ft and 3300ft films for the exact sum curves. The value 
of T c q — 3.6K quoted in the paper was used, and a value of H C 2(0) = lOOfcOe was estimated from the experimental 
plots. There is a clear upturn in H C 2(T) at low T, but it is not explained by the MEF curves which are effectively all 
upturn at these parameter values. Again the fit to exact sum curves seems better. 

In conclusion it seems that the experimental curves fit better to the exact sum curves than the MEF curves. Any 
low temperature upturns in H C 2(T) are not explained by the MEF formula, and presumably correspond to different 
physics. 




FIG. 6. Fitting of Graybeal and Beasley's datau to the MEF formula (solid lines) and perturbation theory with exact sums 
(dashed lines). In each case only a single fitting parameter is needed for all curves, corresponding to the initial slope of the 
TciRa) curve. The films have resistances « 0, 110, 400, 600a 
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FIG. 7. Fitting of Okuma et al's datalij to the MEF formula (solid lines) and perturbation theory with exact sums (dashed 
lines). Different values of the fitting parameter are needed for the two films to get the correct value of T c for the given Ra. The 
films have resistances Ra = 400, 600S1. 
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0.5 1.0 
T/T 

FIG. 8. Fitting of Hebard and Paalanen's dataH to the MEF formula (solid lines) and perturbation theory with exact sums 
(dashed lines). Different values of the fitting parameter are needed for the three films to get the correct value of T c for the 
given Eu. The films have resistances Ea = 2250, 2900, 3300f2. 



10 



ACKNOWLEDGEMENTS 

We thank I. Alcincr, A. Clerk, A.M. Finkel'stein, I.V. Lcrner and Y. Oreg for helpful discussions. RAS ac- 
knowledges the support of a Nuffield Foundation Award to Newly Appointed Lecturers in Science and Mathematics. 
BSH acknowledges the support of a UK EPSRC Graduate Studentship. VA is supported by the US National Science 
Foundation under grant DMR-9805613. 



1 S. Maekawa and H. Fukuyama, J. Phys. Soc. Jpn 51, 1380 (1982). 

2 H. Takagi and Y. Kuroda, Solid State Commun. 41, 683 (1982). 

3 A.M. Finkel'stein, JETP Lett. 45, 46 (1987). 

4 R.A. Smith, MY. Reizer and J.W. Wilkins, Phys. Rev. B 51, 6470 (1995) 

5 A useful review of the whole area can be found in: A.M. Finkel'stein, Physica 197B, 636 (1994). 
6 H.R. Raffy, R.B. Laibowitz, P. Chaudhari and S. Maekawa, Phys. Rev. B 26, 6607 (1983). 

7 J.M. Graybeal and M.R. Beasley, Phys. Rev. B 29, 4167 (1984). 

8 D.B. Haviland, Y. Liu and A.M. Goldman, Phys. Rev. Lett. 62, 2180 (1989). 

9 S.J. Lee and J.B. Ketterson, Phys. Rev. Lett. 64, 3078 (1990) 

10 P. Xiong, A.V. Herzog and R.C. Dynes, Phys. Rev. Lett. 78, 927 (1997). 

11 J.M. Valles Jnr, R.C. Dynes and JP. Garno, Phys. Rev. B 40, 6680 (1989); Phys. Rev. Lett. 69, 3567 (1992). 

12 J.M. Valles Jnr, S-Y Hsu, R.C. Dynes and J. P. Garno, Physica 197B, 522 (1994). 

13 J.A. Chervenak and J.M. Valles Jnr, Phys. Rev. B 51, 11977 (1995). 

14 A.F. Hebard and M.A. Paalanen, Phys. Rev. B 30, 4063 (1984). 

15 S. Okuma, F. Komori, Y. Ootuka and S. Kobayashi, J. Phys. Soc. Jpn. 52, 3269 (1983). 

16 R.A. Smith and V. Ambegaokar, in preparation. 

17 S. Maekawa, H. Ebisawa and H. Fukuyama, J. Phys. Soc. Jpn. 52, 1352 (1983). 

18 A. A. Golubov and V.V. Dorin, J. Low. Temp. Phys. 78, 378 (1990). 

19 B. Brandow, Phys. Rep. 296, 1 (1998). 

20 A.M. Finkel'stein, Sov. Phys. JETP 57, 97 (1983); Z. Phys. B 56, 189 (1984). 

21 Y. Oreg and A.M. Finkel'stein, Phys. Rev. Lett. 83, 191 (1999). 

22 A. A. Abrikosov and LP. Gor'kov, Sov. Phys. JETP 12, 1243 (1961). 

23 S. Hikami, Phys. Rev. B 24, 2671 (1981). 



11 



